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Abstract 

We introduce a path-integral approach that allows to compute 
charge density oscillations in a Luttinger liquid with impurities. We 
obtain an explicit expression for the envelope of Friedel oscillations in 
the presence of arbitrary electron-electron potentials. As examples, in 
order to illustrate the procedure, we show how to use our formula for 
contact and Coulomb potentials. 
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In the last years there has been much activity addressed to the study of 
condensed matter and statistical mechanics problems through field-theoretical 
methods In particular the physics of one- dimensional (Id) systems of 
strongly correlated particles has become a very interesting subject since one 
can take advantage of the simplicity of the models at hand and, at the same 
time, expect to make contact with experiments. For instance, the recently 
built quantum wire || is a good realization of a Id electron gas. From the 
theoretical side, the simplest formulation of a Id electronic system is given 
by the Tomonaga-Luttinger (TL) model || which has been successful in de- 
scribing some qualitative features of a Luttinger liquid such as spin-charge 
separation and non-universal exponents in the decay law of correlation func- 
tions |Q. There are, however, two crucial issues that are not considered in 
the original versions of this model: the presence of a non trivial interac- 
tion between electrons and impurities []5| and the effect of long-range (LR) 
electron-electron interactions || p]. As it is well known, the former leads 
to the occurrence of Friedel oscillations in the charge density profile, at least 
for Fermi liquids ||. On the other hand, as the dimensionality of a system 
decreases, charge screening effects become less important and the LR inter- 
action between electrons is expected to play a central role in determining the 
properties of the system. In fact, from a theoretical point of view, the ef- 
fects of LR interactions have been recently discussed in connection to several 
problems such as the Fermi-edge singularity ||, the insulator- metal tran- 
sition and the role of the lattice through umklapp scattering and size 



dependent effects [llj. Thus, it is quite interesting to study the interplay 
between impurities and LR interactions by considering Friedel oscillations 
in a Id system. Some time ago, Egger and Grabert |12| analyzed this phe- 
nomenon. By combining the techniques of standard bosonization jn| with 
the self-consistent harmonic approximation [[14] and quantum Monte Carlo 
simulations they were able to get explicit results for both weak and 
strong impurity scattering regimes. Later on, the authors of Ref. |TB| used 
bosonization and a scattering description to get some exact results for the 
short-range case and for a special value of the coupling constant, equivalent 
to the so called " Toulouse point" in the anisotropic Kondo problem. More re- 
cently, the authors of Ref. [|17] used again standard bosonization to address 
the same problem emphasizing the equivalence between the TL model in 
the presence of a single non-magnetic impurity and a boundary Sine-Gordon 
model. 
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In this work we present an alternative approach to this problem based on 
a path-integral bosonization technique previously developed in the context 
of non local quantum field theories [|18|| . This method seems to be specially 
adequate to consider LR interactions. Indeed, it has recently provided a 
straightforward derivation of the electronic Green's function in the presence 
of non contact potentials |19fl . Then, our main purpose here is to show how 
to extend this formulation to the computation of Friedel oscillations. 



We start from a modified non local Thirring model |20[ described by the 
following (Euclidean) Lagrangian density: 

C = tV(0 + lo k F )y + J d 2 yJ^x)U w {x,y)J^y) + ^ - M(x)^ (1) 

where x = (r x ,x) = (xo,x\), and J M = ^7^. The functions U^(x,y) are 
forward-scattering potentials. Setting C/( ) — ^(1) — S 2 (x — y) one gets 
the covariant and local version of the Thirring model usually studied in the 
context of (1+1) QFT's. 

On the other hand, the choice (x, y) = C/(|x— y|) 5(r x — r y ) and Um(x, y) = 
0, yields the simplest version of the Tomonaga-Luttinger (TL) model with 
an instantaneous distance dependent potential and no current-current fluc- 
tuations. The last two terms in (HD correspond to forward and backward 
electron-impurity scattering, respectively. 

The main purpose of the present paper is to evaluate the v.e.v. of the charge- 
density: 

(p(x)) = + e~ 2ikFX ^ L + e 2 *^^), (2) 

for an arbitrary electron-electron potential Ui^x, y). Also we will be spe- 
cially interested in case that the impurity terms are Cq(x) = V5(x — d) = 
M(x) and C\(x) = 0, where V is a constant proportional to the impurity 
tunneling barrier situated at x = d. Using a suitable representation of the 
functional delta and introducing an auxiliary vector field (see Ref. 
for details), the partition function of the model under consideration reads 



Z = N JVA^ e- s ^det (0 + V2fi + + 4 + i lo k F - M{x) + s^ v e 2lkF ^ 

(3) 
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where S[A] is the free quadratic action for A^, To = I and Ti = 75. Equation 
(H) can be obtain by functional derivation of Eq. (|[) with respect to the 
source s . 

As it is known, the massive-like determinant in equation (0) cannot be exactly 
solved, even in the local case. However, we can take advantage of the fact 
that the vacuum to vacuum functional can be written in such a way that 
non local terms are not present in the determinant. Therefore the terms 

can be decoupled from fermions by performing chiral 
and gauge transformations in the fermionic path-integral measure. Indeed, 
decomposing A^ix) in longitudinal and transverse pieces 

A„(x) = e^d v - + d^{x) - -^(^(x) + s^x)), (4) 

where $ and rj are boson fields (to be associated to the normal modes of the 
system) and applying, as anticipated, functional bosonization techniques |L8 



to express the fermionic determinant in terms of $ and 77, one finally obtains 
Z = Nf VxD X V$Vr] exp (S bos ) exp {-S fer ) exp {-S[M, s„]), (5) 
where Sf er corresponds to free massless fermions (x and x) an d 

S[M, s^} = J d 2 xx (s^x) e 2lk ^ - M(xj) e~ 2 ^ X - (6) 
Concerning S& os , it can be more briefly described in momentum space: 

H-p) 
fj(-p) 

cu-p) 



s bos = I A (Hp) m C'(p))\ ^ B(p) ^ II v(- P ) | + 



(2tt) 2 V = ^fj, 



A(p) 


Cjp) 


E(p) 
2 


C(P) 

2 


B(p) 


EM 


E(p) 
2 


EM 
2 


D{p) 



(27r) 2ly (0) W w (0V 

(7) 

where we have defined C' = + and A(p), B(p),C(p), D(p), E{p) and 
F(p) are potential dependent functions (See J20| for more details). 
At this point we see that the generating functional can be formally expanded 
in powers of — s M e^e 2lfcFxr,/ ^ , in complete analogy with the usual proce- 
dure employed in the path-integral bosonization of (1 + 1) massive QFT's |H 



In fact, the x dependence of this perturbative parameter, together with 
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the appearance of C (x) in the bosonic action are two of the new features of 
the present computation. As far as these functions are well-behaved one can 
assume the existence of every term in the corresponding series. 
From now on we will specialize the computation to the case sj = 0. This 
allows to define M±(x) = M(x) — SQ(x)e T2lkFX and one can then perform 
the above mentioned expansion of Z taking M±(x) as perturbative param- 
eters. As explained in Ref. |2(J one can show that the same expansion can 
be obtained by starting from a purely bosonic non local extension of the 
sine-Gordon model given by 

1 lr 

£' = oO^V?) 2 + o / d 2 yd^(x) d M (x,y) d^{y) + 



2 V pr/ 2 

+ F » ~ ^2 (*+•(*) elMx) + «-(*) e ~ lMx) ) (8) 

where F^(x) represents a couple of classical functions to be related to the 
C'^s and d^{x,y) are two bilocal functions that will be associated to the 
electron-electron potentials (a similar non locality in the kinetic term was 
considered in the study of the influence of LR correlations in the metal- 
insulator transition fl23j). f3 is a constant and a±(x) are functions that can be 



considered as extensions of the parameter ao used by Coleman |24| . Indeed, 
for dfjj) = = Fp and a + = ct_ = «o = constant, the model above coincides 
with the usual sine-Gordon model. In the present approach the quantities 
a±(x) are related to M±(x), which are in turn connected to the strength of 
the scatterer. Let us stress that in our formulation it is straightforward to 
consider a non point-like impurity (a±(x) ^ V5(x — d)). However, in order 
to illustrate our method, in this paper we will consider the usual case of a 
completely localized impurity. For this particular case Eq. (||) contains the 
same terms that can be derived from standard bosonization (See for instance 

[0)- 

Now, going to momentum space and employing standard procedures to 
evaluate each v.e.v., the partition function Z' corresponding to this general- 
ized sine-Gordon model coincides with Z provided that the following three 
relations hold: 



^ (9) 



I bo^(i) (p) + Pi */(<)) (p)] + P 2 2 [p 2 + d {0) (p)pg + d {1) (p)p 2 } 
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= M± (x) = M{x) - s (x)e^ kFX (10) 
2iC' li (-p)e lu ,p v _ pF^-pjp^ ^ 



I \PqU {1) (p) + p\ U {0) (p)] + p 2 p 2 + d (0) p 2 + d (1) pf 

Therefore, we have obtained an equivalence between the partition functions Z 
and Z' corresponding to the non local Thirring and sine-Gordon models with 
extra interactions defined above. This means that we can use Z' together with 
the above conditions in order to compute the charge-density in the Luttinger 
liquid in the presence of impurities. Indeed, as a result of this bosonization 
technique we can evaluate < p(x) > through functional derivation of Z' 
instead of Z. In so doing we obtain 

(p{x)) = (-^=d x ip(x) + cos (2^if(x) - 2k F x)) (12) 



where the v.e.v. is taken with respect to the Lagrangian density C'[so = 0] 
obtained from (|8|) after using equations (|), fllPl) and ([TT|) and setting s = 0. 
Note that we have also set (3 = 2^/7r. 

Let us remark that there is an additional contribution to Eq. (O), coming 
from the functional derivative of the normalization constant N'[C'] with re- 
spect to s^. Since this quantity is a constant, its only effect is to shift the 
background value of the charge density. For this reason we have just disre- 
garded it. 

Now we return to our main goal, that is to use the path-integral framework 
depicted above in order to obtain an explicit formula for the charge- density in 
a Luttinger liquid with an arbitrary electron-electron and electron-impurity 
interactions. When one imposes these conditions in C'[sq = 0] one gets a 
Lagrangian density which has an undefined parity as functional of </?. How- 
ever it is much simpler to work with an even Lagrangian since in this case 
all v.e.v. 's of odd functions of ip will vanish. It is easy to see that the trans- 
lation (p(x) — > <p(x) + f(x) yields an even Lagrangian C' even provided that 
the classical function f(x) is -^-independent and its gradient satisfies: 

d x f(x) + - ( dy U{x - y) d y f(r x , y ) + ^-V 5(x - d) = 0. (13) 

71 J " \/TT 
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We then get 



(p{x)) = -^<9 X / + cos(Vi^/(x) - 2A; F x)(cos(v / 4^(x))> £ ,_, (14) 

'7T 



where cos(v47r/(x) — 2/cpx) is called the Friedel oscillation and A(x) = 
{cos(\/4ft(p(x)))ci even is the corresponding envelope. 

Let us point out that Eq. ( O ) has been previously found in Ref. [T7|]. As 
shown by these authors, in the short-range case it has the solution / = 
constant oc U whose only effect is to add a constant phase in the cosine term 
associated to the Friedel oscillation. A non trivial phenomenon takes place 
for LR potentials, since the cosine ceases to be a periodic function. This 
nonperiodicity effect, although weak at large distances, could eventually be 
observed in carbon nanotubes ||25|1 . 



From now on we shall focus our attention on the computation of the envelope 
of the oscillation. Since C even is not exactly solvable, we shall employ the 
well-known self-consistent harmonic approximation [T4\, which amounts to 
replacing C' even by 

^ l/o N2 1 f r n I \tti \o / \ TTlCV) <5(x — d) 2 

£scha = ■zKOpip) +- / rfy^(r I ,x)[/(x-y)9 y ^(r I ,y)+ (p 

(15) 

where m(V) is a constant, related to the impurity strength, to be variation- 
ally determined. The precise relationship between m(V) and V was obtained 



in Ref. |n|. For instance, in the strong-scattering limit, when V is much 
larger than a certain bandwidth, one has simply m = V (See also Ref. ||17|| ). 
Let us now consider the computation of A(x) using this approximation. Per- 
forming a translation in the field <p(x) —>■ <p(x) + a(x), with a(x) a classical 
function, we find A(x) = exp iy / 7ra(x). 

Going to momentum space we see that the Fourier transform of a(x) satisfies 
an integral equation whose solution is 

aM = -jT^«-Wi-?£2M'l ("J 



with 



oc 




I(p ,r)=l d qi (17) 
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where we defined r = |x — d\. 

The envelope of the Friedel oscillation then reads 

Mr) = exp -\ £ in (iM - .'^tl'o) ) • < 18 > 

which is our main formal result. Indeed, formulae ( |17D and (|18|) give an 
analytical expression (exact within the gaussian approximation) for A(r) as 
functional of both the electron-electron potential and the variational param- 
eter m(V). Since the self-consistent harmonic approximation seems to fail 



for weak impurity strength, due to the neglect of interwell tunneling [|T2 
we restrict our analysis to the strong impurity regime. We also consider a 
large distance approximation of Eq. (|T7|) which consists of inserting 1/r as 
infrared cutoff. Let us call I r (po,r) the integral fll?] ) regulated in this way. 
We will examine, as examples, two specific short range and Coulomb poten- 
tials. This, in turn will allow us to illustrate how to use our general formula 
(|T8D for other cases. Moreover, since these problems were previously consid- 



ered in Refs. JT^| and []T7|| by using standard (operational) bosonization, our 



computation will give an independent confirmation by means of a different 
approach. First of all we note that it is convenient to split out the two terms 
of the exponential factor on Eq. (^), such that A(r) = exp (T + W) (r). 
For the simple contact potential U(qi) = U = constant, we get 

T(r) = ln(Ar) -9 (19) 

and 

W(r) = g exp(mg 2 r) Ei(-mg 2 r) - g exp(2g) Ei(-2g) (20) 

where A is an ultraviolet cutoff. We have also introduced the interaction 
constant g = (1 + — Taking into account the asymptotic behaviour of 
the exponential integral function £7, for mg 2 r » 1, one obtains 

A(r) = C(g, A) (2gr)-° exp( — ) (21) 

mgr 

which coincides with Refs. |T2J and [[T7| under the same regime. 
In the Coulombian case, one has U{\x\) — Uj \J\~x\ 2 + b 2 , whose Fourier trans- 
form is U(qx) = 2UKo(bqi), where b plays the role of a lattice spacing. 
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Inserting this expression in T and W, and considering the same regime as 
before we find that W vanishes and 



T(r) 



7T 

2U 



4C7 6 
— In — 

7r 2r 



7T 2 



(22) 



which yields 



,-l(r) =C'(. y -^A)oxp(-^ln J ). 



(23) 



Again, this behaviour is equal to the one previously found in Refs. [E2| and 



In summary, we have described an alternative bosonization approach to 
the computation of charge- density fluctuations. This technique, previously 
originated in the context of QFT's, parallels, in the path-integral framework, 
the operational schemes usually employed in condensed matter applications. 
In particular, we have computed the envelope of Friedel oscillations in a sim- 
ple version of the TL model with a non-magnetic impurity. By combining 
that bosonization procedure and the self-consistent harmonic approximation 
we were able to express the envelope of the oscillations as a functional of 
the electron-electron interaction (see Eqs. ([[7]) and (18)). Finally, as a con- 
sistency check of this formal result, and in order to illustrate our method, 
we considered the long-distance regime for contact interactions and Coulomb 
potentials. Our results are in agreement with Refs. [12 and [17 
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